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Abstract: The measurement of Higgs self-coupling is one of the most crucial physics 
goals at the future colliders. At the LHC, the di-Higgs production is a main way to measure 
the Higgs trilinear coupling. As a complementary to the di-Higgs production, tthh process 
may open a new avenue to measure di-Higgs physics at the LHC and a future 100 TeV pp 
collider or a high energy e'’'e“ collider since the extra tt in the final states may efficiently 
suppress the backgrounds. However, such a kind of process is also controlled by the top- 
Higgs coupling. In this work, we investigate the impact of CP-violating top-Higgs coupling 
on tthh production at the LHC, e'’'e“ and a 100 TeV hadron collider under the current 
Higgs data. Within 2a Higgs data allowed parameter region, we find that the cross section 
of tthh at the LHC-14 TeV, TeV and VHE-LHC/SPPC-100 TeV can be enhanced 

up to 2.1 times the SM predictions. The future precise measurement of Higgs coupling will 
reveal the nature of top-Higgs interaction and improve the sensitivity of the determination 
of Higgs self-coupling through tthh production. 
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1 Introduction 


In 2012, a bosonic resonance with a mass around 125 GeV was independently observed by 
the ATLAS and CMS collaborations at the LHC [ 1 , 2 ]. Up to now, most measurements 
of its properties are consistent with the predictions of the Standard Model (SM) Higgs 
boson [3,4]. However, the Higgs couplings with top quarks and with itself are still vacant 
and need to be verified at the future colliders. 

The tree-level Higgs potential of the SM is given by. 


V = + A((/)V)^ 



( 0 , V + h)'^ , 


which yields the trilinear and quartic Higgs self-couplings. 


^3h = 


3ml 





( 1 . 1 ) 


( 1 . 2 ) 


where the Higgs boson mass ruh — 125 GeV have been well measured in the experiments and 
the vacuum expectation value of the Higgs field, v = ~ 246 GeV. Hence, the 

determination of the Higgs self-coupling can directly test the relation Eq.(1.2) obtained 
from the minimization of SM Higgs potential Eq.(l.l). At the LHC, the main way to 
measure the Higgs self-coupling is the di-Higgs production, which is dominated by the 
gluon-gluon fusion mechanism and has small cross section [5-7]. Among various decay 
channels, the process hh —)• bb'j'j is expected to have the most promising sensitivity due to 
the low backgrounds at the LHC [ 8 ] while the 46 hnal state has the largest fraction [9] . Some 
applications of advanced analysis techniques to di-Higgs production have been proposed 
to improve the sensitivity of bbT~^T~ and bbW^W~ channel [10]. However, one should 
keep in mind that, the main production process of di-Higgs gg —?• hh can also be triggered 
by the top-Higgs coupling itself through the box diagrams. So, any new physics in the 
top-Higgs coupling may significantly affect the measurement of Higgs self-coupling at the 
LHC [11-18]. 
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In the SM, top quark has the strongest coupling to the Higgs boson and is widely 
speculated as a sensitive probe to the new physics beyond the SM. The most sensitive direct 
way of measuring the top-Higgs coupling at the LHC is the associated production of the top 
pair with Higgs boson that has been extensively studied in literature [19-22]. However, 
since the cross section of pp —?■ tth is about 130 fb at 8 TeV LHC [23]), current LHC 
luminosity and analysis are not yet sensitivity enough to observe such a signal. An upper 
limit on the signal strength to be < 3.9 at 95% C.L. limit has been set 

up by ATLAS collaboration through combining h ^ bb and L —)• 77 channels [24]. While 
the CMS collaboration gives a limit of 0.9 < ptih ^ 3.5 by using all search channels [25]. 
The other direct way is to observe the single top associated production with the Higgs 
boson. Such a process has a small cross section (18.28 fb at 8 TeV LHC) but is advocated 
to determine the sign of the top-Higgs coupling at the LHC [26, 27] . Recently, the CMS 
collaboration has presented the result on thj searches in the L —)• 77 channel, and obtained 
a weak bound on the cross section of events with inverted top-Higgs coupling [28]. 

In this work, we will investigate the effect of non-standard top-Higgs coupling in the 
tthh production at the LHC, a e'’'e“ and a 100 TeV hadron collider under the current 
Higgs data constraints. At the LHC, although the cross section of pp —)• tthh is about 
an order smaller than that of pp ^ hh, the additional tt in the final states may suppress 
one order or orders more backgrounds [29,30]. Besides, such a process has a cross section 
monotonically increasing with respect to Higgs self-coupling [ 6 ], which may complement 
the process pp —)■ hh in measuring Higgs self-coupling, particularly for A S> A^m- So, the 
process pp —)• tthh opens a new avenue to measure di-Higgs physics at HL-LHC and a 
future 100 TeV pp-collider. On the other hand, given the limited precision of the LHC, 
an e^e“ collider is crucial to scrutinize the detailed properties of the Higgs boson that 
might uncover the new physics beyond the SM [31]. In e^e“ collision, the double Higgs 
strahlung e'’'e“ —?■ Zhh is considered as the main production process to measure the Higgs 
self-coupling. Besides, it should be noted that the di-Higgs bosons can also be radiated 
off from top quarks through the process e"’'e“ —>■ tthh, or 77 —)• tthh at photon collider, 
where the energetic electron beam to a photon beam is converted through the backward 
Compton scattering [33]. Similar to pp —)• tthh, the processes e'^e~ —)• tthh not only 

involve the Higgs self-coupling but also are governed by the top-Higgs coupling. So the 
investigation of tthh production may lead us to obtain the nontrivial information on the 
Higgs potential and test the top-Higgs coupling [34,35]. 

The structure of this paper is organized as follows. In Section 2, we will briefly intro¬ 
duce the non-standard top-Higgs interaction and set up the calculations. In Section 3, we 
present the numerical results and discuss the effects of non-standard top-Higgs coupling 
in the di-Higgs production at the LHC, e"’“e“ and a 100 TeV hadron collider. Finally, we 
draw our conclusions in Section 4. 
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2 CP-violating Top-Higgs Couplings 

In this study, we parameterize the top-Higgs couplings using the phenomenological La- 
grangian: 

£ = —-^t(cos 6* + sin 0)t/i. (2.1) 

v2 

In the SM, yt takes the value = \/2mtjv and sin0 = 0, with v ~ 246 GeV being 
the vacuum expectation value of the Higgs held. A pure pseudo-scalar interaction can be 
obtained by setting cos 9 = 0. A CP violating interaction is realized if both cos 9 ^ 0 and 
sin0 7 ^ 0. The exact values of these coefficients depend on the specihc model. Here we 
are interested in a model-independent approach to determine the impact of general top- 
Higgs coupling on tihh production. In the following calculations, we dehne two reduced 
couplings: ct = yt cos 9/ytgj^f and q = yt sin 9/ytg;^,^ to discuss our results. 

Since the CP-violating top-Higgs interaction can sizably alter the production rate of 
gg ^ h and decay width of h —)• 77 through the loop effect, the most relevant indirect 
constraint on the couplings ct and q should be from the Higgs data. The signal strength 
of one specihc analysis from a single Higgs boson is given by 

fi — ^ ^ ( 2 . 2 ) 

i 


where the sum runs over all channels used in the experimental analysis. Each channel is 
characterized by one specihc production and decay mode. The individual channel signal 
strength can be calculated by 


[fj X BR]- 
WsM X BRsM]i' 

(2.3) 

and the SM channel weight is 


_ _ Q [crsM X BRsM]i 

Cj WsM X BRsujj 

(2.4) 


where e* is the relative experimental efficiencies for each channel. But these are rarely 
quoted in experimental publications. In this case, all channels considered in the analysis 
are treated equally, i.e. e* = 1. The reduced effective coupling Chgg and Ch-f-y can be 
parameterized through the ct and ct as following [36], 

Chgg - Ct + 2.6c? -h O.llctict - 1) , 

~ (1.28 - 0.28ci)^ + {OASctf . (2.5) 

We confront the effective coupling Chgg and Ch^y with the Higgs data by calculating the 
of the Higgs sector with the public package HiggsSignals-1.3.2 [37]. We choose the 
mass-centered method in the package HiggsSignals. Although the CP-violating top- 
Higgs interaction can contribute to the electric dipole moment (EDM), the bounds on the 
coupling Ct strongly rely on the assumption of Higgs couplings to other light fermions [38] . 
But such light quark Yukawa couplings are generally unobservable at the LHC, we do 
not impose EDM constraints in the study. Other low-energy physics constraints, such as 
Bs — Bg and B —)• are still too weak due to the large uncertainty [36]. 


- 3 - 




3 Numerical Resutls and Discussions 


In our numerical calculations, we take the input parameters of the SM as [39] 

mt = 173.07 GeV, mw = 80.385 , mz = 91.19 GeV, 
mh = 125.9 GeV, sin^ = 0.2228, a{mz)~^ = 127.918. (3.1) 

For the strong coupling constant as{fJ-), we use its 2-loop evolution with QGD parameter 
= 226 MeV and get Ugimz) = 0.118. We use GTIO parton distribution functions 
(PDF) for the calculation oipp —)• tthh [40]. The renormalization scale pR and factorization 
scale pf are chosen to be pR = Pf = {^h + nT-t)- It should be noted that the higher 
order corrections are usually needed to improve the reliability of the leading order results. 
However, such calculations for four-body production involve complicated techniques and 
are beyond our study in this work. In the following we will use the ratios of the leading 
order cross sections to present our results, which have the weak dependence on 

the variation of the scale. Besides, it is expected that the high order effect on the cross 
sections can be largely canceled in the ratios. We perform the numerical calculations by 
the package calchep-3.4 [41]. 

3.1 LHC and VHE-LHC/SPPC 



Figure 1: Feynman diagrams for the partonic process gg —)• tthh at hadron collider. 
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At hadron collider, the main contributions to tthh production are from the gluon fusion 
processes. The corresponding Feynman diagrams for the partonic process gg —)• tthh are 
shown in Fig. 1. The qq annihilation processes can be obtained by replacing the initial 
gluons with qq in the s-channel in Fig. 1. According to the production of di-Higgs bosons, 
we can classify amplitudes of pp —>• tthh into two categories: one is proportional to agy^] 
the other one is proportional to asUt^sh- In the SM, both amplitudes have the same sign 
and are constructive, which lead to the cross section of pp —)• tthh monotonically increases 
with respect to Higgs self-coupling A 3 / 1 . However, it should be noted that each amplitude is 
modulated by the top-Higgs coupling. So the measurement of Higgs self-coupling through 
pp tthh strongly depends on the determination of top-Higgs coupling. 
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Figure 2: Ratios of in the plane of q — ct at 14 TeV LHC and VHE- 

LHC/SPPC, where the dashed contours correspond to the 68% C.L. and 95% C.L. limits 
given by the current Higgs data fitting. 

In Fig. 2, we respectively show the ratios of in the plane of q — ct at 

14 TeV LHC and VHE-LHC/SPPC, where the dashed contours correspond to the 68% C.L. 
and 95% C.L. limits given by the current Higgs data fitting. From Fig. 2, we can see that 
the negative scalar component ct in Eq. 2.1 is strongly disfavoured because the experimental 
measurement of the Higgs diphoton rate are consistent with the SM prediction. The 
pseudoscalar component in the range |q| > 0.6 have been excluded at 95% C.L by the 
Higgs data fitting. Besides, when c* = 1 (or yt = a narrow region of pseudoscalar 

component |q| < 0.3(0.44) is still allowed at 68% (95%) C.L.. 

As mentioned before, the cross section of pp —)• tthh not only depends on the magnitude 
of the top-Higgs coupling but also on the relative phase angle between between ct and c*. 
(i) When q = — 1 and Q = 0 (or yt = —yt^), the ration of is smaller 

than 1 due to the deconstructive interference between the processes with coupling ccsT/f 
and those with coupling asyt^sh] (h) when q = 0 and Q = 1 (or yt = iyf^), there is no 
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interference between the terms with asUt those with asVt^zh: which will reduce the 
cross section of pp —)• tihh. Finally, we find that the value of the ratio 
can maximally reach about 1.5 (2.1) in the 68% (95%) C.L. allowed region at 14 TeV LHC 
and VHE-LHC/SPPC. 




Figure 3: The parton-level transverse momentum distributions of top quark and Higgs 
boson at 14 TeV LHC. 

We implement the CP-violating interaction of tih in Eq. 2.1 with the package FeynRules 
[42]. In Fig. 3, we show the parton-level transverse momentum distributions of top quark 
and Higgs boson in the presence of the CP-violating interactions by using the package 
MadGraphS [43] at 14 TeV LHC. From Fig. 3, we can see that the shapes of pT distributions 
of top quark and Higgs boson in the SM are almost the same as those in the CP-violating 
interactions. 




Figure 4: The parton-level distributions of polar angle AOth and the azimuthal angle A(pth 
between the top quark and the Higgs boson at 14 TeV LHC. 
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Since the CP-violating tth interaction in Eq. 2.1 may affect the polarization states 
of the top quark, we present the parton-level distributions of polar angle A6th and the 
azimuthal angle between the top quark and the Higgs boson. From Fig. 4, we can 

find that the CP-violating interactions slightly change the distributions of and Acpt^, 
as the comparison with the SM. To be specific, when ^ 7r/2 {AOth > 7112), the SM 
has more (less) events than the CP-violating interactions. Similar results can be observed 
in the distribution of A(j)th- Here it should be mentioned that such differences will be more 
obvious if the pseudo-scalar component ct becomes large. However, as shown in Fig. 2, the 
current Higgs data has already tightly constrained the size of q. Besides, these differences 
will be further diluted when the detector effects are taken into account. 

The LHC observability of tihh in the SM including the parton shower and detector 
effects has been investigated through both Higgs bosons decaying /i —)• 66 and semi-leptonic 
and hadronic top decays in Ref. [29]. The main backgrounds are tihbb, ttbbbb and tiZbb 
productions. Due to the four 6-jets in the signal tthh coming from two Higgs bosons, the 
most efficient cut is the reconstruction of the Higgs mass through selecting the Higgs decay 

where i ^ j ^ k ^ I run over 


■ , 1 .... 2 {mi, 

jets by minimizing, Xhh ~ -- 


all 6-tagged jets and A^ = 20 GeV. To remove the kP-|-jets background, at least one top 
quark is required to be reconstructed by using Xt, — ^ - — for semi-leptonic and 


,, {mj. j, ju~mt)^ 

Xti,^ = —- for hadronic top decays, respectively. Since the Higgs data has tightly 

constrained the CP-violating tth interactions (see Fig. 2), the allowed values of ct and ct 
can not sizably affects the polarization states of the top quark and will not change the 
distributions of the final states. Therefore, we can expect that the cuts efficiency obtained 
in the Ref. [29] is applicable to our case, in particular for pure scalar tth interaction. In 
this case, we can estimate the LHC sensitivity of our signals by normalizing the results in 
Ref. [29]. For example, when q = 1.12 and q = ±0.23 (at 95% C.L. Higgs data allowed 
region), = 2 and S/B = 30%. The corresponding statistical significance 

can reach 3cr at 14 TeV LHC with the luminosity L = 3000 fb“^. So we may expect that 
the enhancement of tihh arising from the CP-violating tth interactions may be observed 
at the HL-LHC. 


3.2 e'''e and 77 collisions 

At a high energy e'’'e“ collider, e"*“e“ —)• tihh process can occur only through the s-channel. 
The corresponding Feynman diagrams for e+e“ —)• tihh are shown in Fig. 5. Different from 
the process pp tihh, the di-Higgs bosons can be produced not only through the top- 
Higgs or Higgs self interaction but also by the Higgs gauge couplings hZZ or hhZZ. Since 
we assume that the Higgs gauge couplings and Higgs self-couplings be the SM values, any 
modifications in tih coupling can change the interference behavior between the amplitudes 
with asVt and with asVthh- 

In Fig.6, we show the ratios of in the plane of ct — q at a e+e“ 

collider with ^/s = 1 TeV. From Fig.6, we can see that the cross section of e'’'e“ —)• tihh 
monotonically increases with ct when q = 0 and can maximally reach 4 times as large as the 
SM prediction. For a given q, the pseudo-scalar component in tih can reduce or enhance 
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Figure 5: Feynman diagrams for the process e'’'e —)• tthh. 



Figure 6: Same as Fig.2, but ratios of at a e+e collider with 

Vs = 1 TeV. 
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the cross section of e'’'e“ —>• tthh. Although some of the subprocesses of e^e“ —?■ tthh that 
can be triggered by the Higgs gauge interactions, the cross section of e"*'e“ —>• tthh is still 
dominated by top-Higgs coupling. Within the 68 % (95%) C.L. Higgs data allowed region, 
we note that the value of the ratio can maximally reach about 1.5 

(2.1) at 1 TeV e+e“ collider. 



25 26 27 28 29 30 


Figure 7: Feynman diagrams for the process 77 —)• tthh. 


At a photon collider, the 77 collisions can be also achieved by the inverse Compton 
scattering of the incident electron- and the laser-beam, the number of events can be cal¬ 
culated by convoluting the cross section of 77 collision with the photon beam luminosity 
distribution: 


N^^—).hh / - /lg+g - (T..y.y (S) 


(3.2) 


77 


where is the photon-beam luminosity distribution. a.y.y^hh{s) ( s is the squared 

center-of-mass energy of e'’'e“ collision) is defined as the effective cross section of 77 —>■ hh, 
which can be written as [44] 


/ ^max r^max 

z .s') I F.yyg(x)Fyyg( ) (3-3) 

/a J z‘^/^'max ^ ^ 


Here F^/g denotes the energy spectrum of the back-scattered photon for the unpolarized 
initial electron and laser photon beams, which is given by 


Fy/g(x) 


1 

m 


1 — X + 


1 


1 — X 


4x 4x^ 

^(1 — x ) ^ 2(1 — x )2 


(3.4) 
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with 


1 


m 



~^) + 0 + 2 


8 


2(1 + 0 "' 


(3.5) 


Here ^ = 4£'e£'o/^e (-®e is the incident electron energy, where Eq is the initial laser 
photon energy) and x = E/Eq with E being the energy of the scattered photon moving 
along the initial electron direction. In our calculations, we take the parameters as ^ = 4.8, 
11(0 = 1-83 and Xmax = 0.83 [44]. At photon collider, 77 —)■ tthh process can occur only 
through the t-channel. The corresponding Feynman diagrams for 77 —)• tthh are shown in 
Fig.7. 



Figure 8: Same as Fig. 2, but ratios of at a photon collider with y/s = 1 

TeV. 


In Fig. 8 , we present the ratios of in the plane of q — ct at photon 

collider with ^/s = 1 TeV. Similar to the process pp —?■ tthh, the cross section of 77 —)• tthh 
is determined by top-Higgs coupling. From Fig. 8 , we can see that the allowed maximal 
cross section of 77 —?• tthh can be 2 times as large as the SM prediction, which occurs at 
Ct = 1.2 and q = 0. The pseudo-scalar component in tih always reduce the production 
rate of 77 —)• tthh. 

Table 1: The SM cross sections of signal tthh and the backgrounds tthbb, ttbbbb, ttZbb at 
a e+e” collider with ^/s = 1 TeV. 


e'’'e 

tthh 

tthbb 

ttbbbb 

ttZbb 

0 - (10-3 fb) 

13.55 

6.69 

5.52 

16.24 


10 


















In Table 1, we present the cross sections of e^e“ —?• tthh and the backgrounds e^e“ —)• 
tthbb, ttbbbb, tiZbb for ^/s = 1 TeV. It can be seen that the cross section of tihh can reach 
13.55 X 10“^ fb in the SM due to the suppression of phase space. Given the enhancement 
effect in tihh from the CP-violating top-Higgs couplings, if a e+e" collider could deliver 
an integrated luminosity up to 10 ab“^, there will be about hundreds of tihh events, 
which may be potentially used to test the anomalous top-Higgs couplings. However, the 
exact detector configurations are not finalized. Therefore, the background estimation and 
signal extraction strategies would be largely dependent on the detector designs and trigger 
conditions [45], which is beyond the scope of our study. Besides, we find that the cross 
section of 77 —)• tihh can only reach 3.74 x 10“^ fb in the SM for ^/s = 1 TeV and is hard 
to be detected at a photon collider. 

4 Conclusion 

As a complementary to the di-Higgs production, tihh process may open a new avenue to 
measure di-Higgs physics at the LHC and a future 100 TeV pp collider or a high energy 
e'''e“ collider. In this paper, we studied the impact of CP-violating top-Higgs coupling 
on tihh production at the LHC, a e'’'e“ and a 100 TeV hadron collider under the current 
Higgs data. Within the allowed parameter region, we find that the cross section of di- 
Higgs production at the LHC-14 TeV, e+e"-! TeV and VHE-LHC/SPPC-100 TeV can 
be enhanced up to about 2 times the SM predictions since a large deviation of top-Higgs 
coupling is still allowed. It is expected that the future precise measurement of Higgs 
coupling will reveal the nature of top-Higgs interaction and improve the sensitivity of the 
determination of Higgs self-coupling through tihh production. 
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